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Abstract 

Within the framework of a generalized coUinear factorization for multiple parton 
scattering in nuclear medium, twist-4 contributions to DIS off a large nucleus can be 
factorized as a convolution of hard parts and two-parton correlation functions. The 
hard parts for the quark scattering in the light-cone gauge correspond to interaction 
with a transverse (physical) gluon in the target, while they are given by the second 
derivative of the cross section with a longitudinal gluon in the covariant gauge. We 
provide a general proof of the equivalence of the hard parts in the light-cone and 
covariant gauge. We further demonstrate the equivalence in calculations of twist-4 
contributions to semi-inclusive cross section of DIS in both lowest order and next 
leading order. Calculations of the nuclear transverse momentum broadening of the 
struck quark in the light-cone and covariant gauge are also discussed. 
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1 Introduction 

In the study of nuclear matter and quark-gluon plasma in high energy lepton- 
nucleus, hadron-nucleus and nucleus-nucleus collisions, hard processes such as 
jet and high transverse momentum hadrons production are useful tools as their 
initial production rates can be calculated within perturbative QCD (pQCD). 
Modification of the final jet or hadron spectra known as jet quenching [T] 
due to further interaction between the energetic partons and the nuclear or 
hot medium can be used to probe the properties of medium such as parton 
correlation or gluon density [2]. Indeed, strong jet quenching has been observed 
both in high-energy lepton-nucleus [3] and nucleus-nucleus collisions [IfS] . 
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Current phenomenological studies of experimental data on jet quenching to 
extract properties of the nuclear matter or hot medium rely on pQCD calcu- 
lations of the parton energy loss or modification of the parton fragmentation 
functions due to gluon radiation induced by multiple scattering during the 
parton propagation. Among many theoretical studies [6H7] . twist expansion 
approach p[9] was applied to medium modification of the fragmentation func- 
tions and parton energy loss in both nuclei in deeply inelastic scattering (DIS) 
and hot QCD matter in heavy-ion collisions. Such an approach is based on the 
generalized factorization framework for multiple parton scattering in nuclear 
medium first developed by Luo, Qiu and Sterman (LQS) [TU]. Within this 
framework, LQS proved that the leading twist-4 contribution from multiple 
parton scattering in a large nucleus can be factorized as a convolution of hard 
parts and two-parton correlation functions. Such parton correlation functions 
have leading contributions that are proportional to the nuclear size Ra ~ A^^^ 
and therefore the twist-4 contributions are enhanced by the nuclear medium 
as compared to the same contribution in nucleon collisions. 

The twist expansion technique has been used to calculate the effect of multiple 
parton scattering in dijet and photon production in DIS off nuclear targets 
PTfT^ and Drell-Yan (DY) dilepton spectra inp + A collisions [TSUI]. In these 
calculations, processes involving secondary scattering with a soft gluon (soft 
scattering) are treated in the framework of coUinear expansion in a covariant 
gauge while processes of secondary scattering with a hard gluon (double hard 
scattering) are calculated directly in the light-cone gauge with coUinear ap- 
proximation (neglecting the transverse momentum of the initial gluons from 
the nucleus). In such separate treatments of soft and double hard scattering, 
it is difficult to include the interference between the two type of processes. For 
large values of the final transverse momentum of the dijets, photons or DY 
dilepton, one nevertheless can neglect the interference effects since the forma- 
tion time of the particle production, which dictates the interference, is much 
smaller than the nuclear size. In the investigation of nuclear modification of 
parton fragmentation functions and induced parton energy loss pi9] , it is crit- 
ical to consider soft gluon bremsstrahlung. Since the formation time of soft 
gluons can be larger or comparable to the nuclear size, one has to include the 
interference between soft and double hard scattering which is equivalent to the 
Landau- Pomeranchuck-Migdal (LPM) interference [13]. In this case, both soft 
and double hard processes and their interference terms should be calculated 
within the general framework of coUinear expansion in the same gauge. 

In a covariant gauge, the longitudinal component of the gauge field ^4+ is 
considered large. One usually makes a coUinear expansion of the hard part 
of quark and (longitudinal) gluon interaction. These hard parts do not corre- 
spond to any physical quark-gluon scattering. If one expands the hard parts in 
the transverse momentum of the longitudinal gluon field, the coUinear terms 
are directly related to the hard parts of the vacuum diagram without final 
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gluon interaction. They only contribute to the gauge hnk of the initial quark 
distribution function of the vacuum diagram. The higher order terms in the 
collinear expansion will contribute to higher twist cross section involving two- 
parton correlation functions and derivatives of the corresponding hard parts. 
In a light-cone gauge A'^ = 0, contributions to the higher twist cross section 
can be written directly as a convolution of the two-parton correlation func- 
tion and the hard part of collinear quark and (transverse or physical) gluon 
interaction. Since the two-parton correlation functions can be expressed in a 
gauge-invariant form, the collinear hard part in the light-cone gauge should 
be equivalent to the derivatives of the hard part in the covariant gauge and 
the final higher-twist cross section should be gauge independent. 

In this paper, we first provide a general proof of the equivalence of calculations 
in both light-cone and covariant gauge of the leading twist-4 contributions to 
multiple parton scattering processes at any order of in DIS off a large 
nucleus. We then demonstrate the equivalence with explicit calculations of 
lowest order semi-inclusive DIS cross section, nuclear transverse momentum 
broadening and induced gluon radiation via secondary parton scattering in 
DIS in both light-cone and covariant gauge. 



2 Generalized collinear factorization 



For simplification, we only consider the semi-inclusive differential cross section 
for the final quark and do not consider quark fragmentation. The differential 
cross section for e{l) + A{p) e{l') + + X can be written as 

nsQ^ d^£± El' 



where Lf'^ilJ') = 4 [/^^r + H'^ -(/■/') c/'^l is the leptonic tensor, p = [p+,0,Oj 
is the longitudinal momentum per target nucleon, q = [—Q^/2q~,q~,0] the 
4-momentum of the virtual photon, s = (/ and £ is the 4-momentum of 

the outgoing quark. 

In a large nucleus or hot QCD matter in heavy-ion collisions, a produced ini- 
tial parton may experience additional scattering with other partons from the 
medium. The additional scattering may further induce additional gluon radia- 
tion and cause the parton to lose energy. Such multiple scattering and induced 
gluon radiation will effectively lead to modification of the parton fragmenta- 
tion functions in a medium, photon and dilepton production. These additional 
contributions from multiple scattering are higher-twist contributions and are 
always power- suppressed. They generally involve high- twist matrix elements 
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Fig. 1. Two-gluon interaction in DIS 



of parton correlation of the medium. For multiple scattering in cold nuclear 
matter, two-parton correlations can involve partons from different nucleons 
in the nucleus, they are proportional to the size of the nucleus and thus are 
enhanced by a nuclear factor A^^^ as compared to two-parton correlations in 
a nucleon. 

In general, the leading (or nuclear enhanced) twist-four contributions at all 
orders in to the hadronic tensor of DIS off a nucleus involving a secondary 
scattering with another gluon (see Fig. [1]) can be written as 



J {2'kY J 27i¥J 27iY J ^ J ^ J ^ 



(27r)4y (27r)4y (27r 



X J-Tr[H^^'^^{k,h,hA{A I mtMyi)My2)m I A) (2) 



where Hj^^lk, ki, is the perturbative hard part of the multiple parton 
scattering at any order of the strong coupling constant, n = [0,1, 0_l] and 
n = [1,0, Ox] are unit four-vectors. Here we have suppressed the color index. 
Summations over color indices of the field operators in the matrix element 
and average over the color indices of the initial state partons in the hard 
part are understood. To simplify the notation, we will sometimes suppress 
the polarization indices associated with the photon interaction when they are 
not necessary, H^^i^k, ki, ^2) = H^l{k, ki, k^)- In the above expression, one can 
always assume collinear approximation for the quark fields, k = xp. The trans- 
verse momentum of the quark will only contribute to high twist contributions 
that are not nuclear enhanced. We can further neglect the minus-components 
of the gluon's momentum, which also contribute to higher-twist terms that 
are not nuclear enhanced. One can therefore express the gluon's 4-momentum 
as ki = xip^ + = + k^^. One can also decompose gluon fields as 

AP = A^_^+A''_^+A''_. Here A'[ = A^df"" , A% = pPA-n/p-n, A^_ = nPA-p/p-n, and 
dpa = —Qpa + fipUfj + Upficr IS the trausverse projection [dp^APB" = A±_ ■ B±). 
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2.1 Light-cone gauge 



Let us first consider the above hadronic tensor in the light-cone gauge , A^. = 0. 
In this case, A'L component contributes only to higher twist (beyond twist-4). 
Therefore, the matrix elements of the non-perturbative part of the hadronic 
tensor only involve one term; 



{A I ijmA'[{yM1iy2my) I A) 

= -U^-^iA I miMvi) ■ My2)^{y) I (3) 

for unpolarized nuclei. 

Taking collinear approximation (neglecting gluons' transverse momentum), 
one can carry out all the integrations except y~ and kf = Xip'^. One can 
further integrate by parts in to convert —Xip'^A±{yi) ■ X2P~^A^{y2) 

into gluon field _strength, d+A±{yi) ■ d+A±{y2) = F+^{y^) ■ F+_l(|/^) (Note 
that = dj^A_[_ in the light-cone gauge, A'^ = 0). The hadronic tensor in 
the light-cone gauge can be reduced to 



w(')= [ f dxdxidx2e'''^''y^^'^'^^y^+'^^^y~ 

X {A I mp-^AvT) ■ F^±{y2 m-) I A) 

X Tr[^'"^(xp, xip, X2P)4^], (4) 

Ax 1X2 



where we have suppressed both the + and transverse components of the co- 
ordinates when they are set to zero, = 0, y± = 0, in variables of the field 
operators. 

The above formula is normally used to calculate contributions from double 
hard scattering where gluons carry finite momentum fraction Xi and X2- How- 
ever, in soft scattering processes, the momentum fractions Xi and X2 go to 
zero. One should then choose a regularization prescription on these light-cone 
poles corresponding to specific boundary conditions for the gluon field [16]. In 
principal, all the prescriptions are equivalent for the complete set of diagrams. 
Alternatively, one can calculate both types of multiple parton scattering and 
their interferences in the covariant gauge as was outlined by LQS [TU] . 
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2.2 Covariant gauge 



In the covariant gauge, both the transverse and plus components of the gluon 
fields contribute to the same twist operators while the minus component 
contributes only to higher twist terms. The matrix elements of the non- 
perturbative part can therefore be expanded into four terms, 



{A I m{KiyT)Aliy^)+Kiy^)Aliy^) + Al(yr)^+(y2~) 

+ Ai{yi)Al{y^)My~)\A). 

To isolate the leading terms of the twist-4 contribution with nuclear enhance- 
ment, one can follow LQS [10] to make the following coUinear expansion of 
the perturbative hard part in terms of the gluon transverse momentum , 



dHP''{k,xip, k2) 
dk^ 

52 ZJPCF 



I, k°' 
\K2=X2P "'2± 



dkfdk^ |fci=a;ip,fc2=x2P t^i±t^2± 



(5) 



where terms d^HP"" /dk'^dkf and d'^HP" /dk^dk^ are not listed as they don't 
contribute to the nuclear enhanced twist-4 terms. The first (coUinear) term can 
be reduced to the eikonal term of the twist-2 contribution making the twist-2 
quark distribution function gauge invariant. The rest of the expansion will 
contribute to higher twist terms. Again, after integration by parts and con- 
sidering d+A^^d+Al = d+A± ■ d+A±dP''/2, d+A'[dlA+ = d+A± ■ d±A+dP''/2, 
dj_A^d'^A^ = d±A^ ■ d±A+dP" /2 in a unpolarized nuclear target, one gets 



1^(2) 



1 fdy dy.dy^ I ^^^^^^^^^ix,p+y,+^X2P^y,+^xpry- 



^aP^^ 



2ti 2ti 2tx 

d'^Hp^{k,ki,k2] 



■j/pPp'^ 



- ki=xip,k2=X2P 



+ d'^^Ti 



dkfdk^ 

x{A I m^d^A+iy^) ■ d^A+{y^)ij{y-) \ A) 
dHp„{k,xip, k2) jjpP 

k2=X2P 



dk? 



Xi 



^ Note that when only the relative transverse momentum k± = ki± 
considered as in Refs. [8], H{ki± — k2±) = H{2k±) 



k2± is 
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x{A I mid+A^iVi) ■ d^A+{y^)ij{y-) \ A) 
dHp„{k, ki,X2p) jjp" 



dkf 



X2 



ki=xip 

x{A I mid^Mvl) ■ ^+My^)^p{y~) \ A) 



X1X2 



x{A I ^(0)7i9+Ax(yr) . d+A^{y:^)i,{y-) \ A) }. 
Using the following identities which we will prove in the next section, 



(6) 



d 



'al3 



Q2 



dkfdk 



Hp„{k, ki, k2)fp'' 



ki=xip,k2=X2P 
1 



= dP^Hp^i^k.Xip, X2p) 

X1X2 

d - 1 
= -d^^'^r^Hp^ik, xip, k2) \k2=x2P P''— 

Ok2 Xi 

d 1 
okf X2 

one can reorganize the hadronic tensor in the covariant gauge as 



(7) 



8 J 



dy dy^ dy2 



X 



27r 27r 27r 

d'^HP"{k,ki,k2 



dxidx2dxe''''P^y^ +i^2P+y2+i^p+y- 



dkfdk^ 



-llPpPa 



{A I 



J ki=xip,k2=X2P 

{d±A^{yi) ■ d^A^{y2) - d+A^iVi) " ^±^+(^2") 
-d±A^{yi) ■ d+My2) + d+My^) ■ d^A^{y^)} ij{y 
dy- dyi dy2 



2tt 27r 27r 



-Tr 



d'^HP'^jk, ki,k2 
dkfdk§ 



-iIPpPct 



ki=xip,k2=X2P 



A) 



x{A I ij{0)LF^^{y-) . I A) , (8) 

where we have neglected higher order corrections to the definition of the gluon 
field strength in the covariant gauge. In the original derivation of the twist 
expanion framework ^Q\, LQS neglected terms that are proportional to the 
transverse gluon fields in the covariant gauge in Eq. ([6]), arguing that these 
contributions are suppressed by As we have shown in the above, the 

corresponding hard partonic parts of all these contributions are equivalent, 
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Fig. 2. Diagrammtic Ward identity (dash line represent the cut) 



allowing us to combine all terms. The final result is proportional to the corre- 
lation of total field strength with one common hard partonic parts. Since the 
gluon field strength correlator is now in a gauge invariant form (with gauge 
links from other soft gluon interaction), the hard partonic parts should also 
be gauge invariant and can be calculated in any gauge. 

Comparing Eqs. (jl]) and ([8]), the leading nuclear enhanced twist-4 contribu- 
tions to the hadronic tensor of DIS off a large nucleus in both light-cone 
and covariant gauge have the same non-perturbative matrix elements for two- 
parton correlation functions. Since these matrix elements are gauge invariant, 
the corresponding hard parts should be equivalent and also gauge invariant, 
therefore can be calculated in any given gauge. The equivalence of the corre- 
sponding hard parts can be proved in all orders as given by the identity in 
Eq. dZD- 



3 Equivalence of hard parts 



In this section we will prove the identities in Eq. ([7]). The technique of Ward 
identity has been used frequently in the proof of factorization in pQCD hard 
processes in which longitudinal gluons can be factorized from the hard part 
and give rise to eikonal lines in hard scattering [T7]. This method is also 
the basic ingredient in our proof of the identities that are used to prove the 
equivalence of the hard parts in the leading twist-4 contributions to DIS off a 
large nucleus in the light-cone and covariant gauge. 

We first consider a right-cut diagram with a single external gluon attached to 
the hard part in Fig. [2l The gluon momentum is = Xip'^ + k^^. We again 
neglect the minus component of the gluon momentum which only contributes 
to higher twist terms. Assuming the total amplitude that includes all 
possible gluon attachment to the hard part, gauge invariance, Xipkl = 0, 
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leads to the following Ward identity, 



Hp{k, ki)k{u{k - ki) = -H{k)-t^iu{k - ki), 



(9) 



as shown diagrammatically in Fig. [21 where u{k) is the quark spinor, Hp{k, ki) 
represents the hard part with one gluon attached after the photon coupling and 
H{k) represents the hard part without any gluon attachment. Using equation 
of motion — }^i)u{k — ki) = 0, we get, 



Hp{k, ki)klu{k - ki) = H{k)u{k - k^) 



(10) 



Truncating the quark spinor which sits inside the non-perturbative parton 
matrix elements, we have the following general identity which is also valid for 
the diagram with the off-shell initial-state quarks. 



Hp{k,ki)kl = H{k) 



'111 



Making a collinear expansion in ki = xip'' + fc^j^ of the left side. 



Hp{k, ki)kl = Hp{k, xip)xipP + Hp{k, Xip)k1^ 

dHp{k,ki) 



+ 



dkf 



k'^^xipP + .... 



:i2) 



fci=xip 



and using Eq. f|TT]) for both ki and Xip, Hp{k, xip)xip'' = H{k) = Hp{k, ki)ki, 
we obtain the following identity. 



dHp{k, ki 



dkf 



kf^xipp = -Hp{k,xip)k1 



(13) 



k\=xip 



Note that authors of the Ref [18] have recently derived the above identity case 
by case. 

For diagrams with two attached gluons that contribute to twist-4 terms, one 
can derive two identities similar to Eq. f|TT]) by contracting with the momentum 
of gluon 1 and gluon 2, ki, /c2 and then make collinear expansion in them 
respectively. 



dHp„{k,ki,k2) 



dkf 



kf^XipP = -Hpa{k, Xip, k2)k{^ 

ki=xip 



(14) 
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dHp^{k, ki,k2) 



k21_X2P'' = -Hpa{k, ki,X2p)k2j_ 

k2=X2P 



(15) 



Making again coUinear expansion in k2± on both sides of Eq. (iMl) . multiplying 
both sides with X2p'^ and using Eq. f|T5l) with ki = Xip, one has 



d^Hp^{k, ki,k2 



dkfdk^ 



ki^k!2j_xip^X2p'' 

ki=xip,k2=X2P 

dHpa{k,xip, fcs) 



dkl 



k^^k1^X2p'' 

k2=X2P 



Hp^{k,xip, X2p)k{^k2^ 



(16) 



Similarly making coUinear expansion in ki_[_ on both sides of Eq. (IT^ . multi- 
plying both sides with xipf and using Eq. (JH]) with k2 = X2P gives, 



d^Hp^{k, ki,k2) 



dkfdk^ 



ki=xip,k2=X2P 



k?, kt 



dHp^{k,xip, k2 



dkf 

Hpa{k, xip, X2p)k^j_k2_i_ 



fci=x'ip 



(17) 



We can further integrate the above two sets of equations over the azimuthal 
angle of the gluons' transverse momentum, divide both sides by X1X2 and 
obtain finally the identities in Eqs.([7]), 



dkfdk 



P PpPa 



dP''Hf"'{k,xip, X2P)- 



ki=xip,k2=X2P 



X1X2 



dHP''{k,xip, k2) Pp 



dk2 xi 
dHP''{k,k^,X2p) p„ 



dkf 



X2 



k2=X2P 



k\=xip 



(18) 
(19) 
(20) 



4 DIS in leading order 



The general proof of the identities for the hard parts in the previous section 
is valid for all orders in a^. In this and next section, we will calculate nuclear 
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modification to the semi-inclusive cross section of DIS off a large nucleus up 
to 0{as) as an explicit example to demonstrate the gauge invariance. 



(a) 



(b) 



(c) 




Ap 



Ap 



R 







^1 I ' 





yj 


.V 





Ap 



R 



1 .1 1 c 
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yi \ 


yi y 



Ap 



Fig. 3. DIS in lowest order with zero, one and two gluon exchange between the final 
quark and the target remnant 

The semi-inclusive cross section off DIS in Eq. ([T]) at the lowest order in 
has been calculated in Ref. [19] up to two gluon exchanges as shown in Fig. [31 
The semi-inclusive hadronic tensor can be expanded in terms of the number 
of physical gluon exchanges [19] Wij_^{q,p, i) = ^IiKq^P, 



dHi 



dHi 



1 
1 



d^k 

d^k d^ki 



S^'\i^-k^)TT[Hj^J{x) ¥'\k)]; 



27rJ (27r)M2vr)4 
xTTmf'>Pix,Xr)u/¥'\k,, 



d ki d k2 d k 



c=L,R,C 



271 J (27r)4 (27r)4 (27r) 
xTr[^(2/)^'^(x, x,,X2)u;f'u;/^^Z'i^u k^, k)] 



(21) 



(22) 



(23) 
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where the summation is over different cut diagrams with given number of gluon 
exchanges, Hj^J, Hj^j^^^ and H^^j)^^^" are the hard parts of the corresponding 

cut diagrams (see Ref. [19] for details). For single gluon exchange, kR_[_ = k±, 
kL± = k± + kij_ For two gluon exchange diagram, kRj_ = k±, kci_ = fc± + 
and kii. = k± + ki_[_ + A;2±. The projection operator is defined such that 
i^pkia = {ki — Xip)p. The above expression for semi-inclusive cross section is 
obtained through coUinear expansion in {ki — Xip) of the partonic hard parts 
and the decomposition of the gauge field Ap = ppA'^/p'^ +ujpAa. Generalized 
Ward identities are used to relate the derivatives of the partonic hard parts to 
the hard parts with extra longitudinal gluon attachments. Consequently, the 
unintegrated parton matrix elements contain contributions from all Feynman 
diagrams with different number of unphysical gluon exchanges between the 
propagating quark and the nucleus. They can be defined in a gauge invariant 
way as. 



mik) = j t/V^-^(^l^(o)/:(o,y)^(y)|A), (24) 
l>W(^i,^) = / AAie^'-^+^'=^-^MA|^(o)/:(o,yi)z},(yi)/:(yi,y)^(y)|A), 

(25) 

^^^}{ki,k2,k) = j (fyd''yi(fy2e'^-y^'^'-^'^'^'-''' 

X (A|^(0)£(0, y2)Dp{y2)C{y2. y,)D,{y{)C{y^, y)i'{y)\A), 

(26) 

where Dp{y) = dp + igAp{y) is the covariant derivative and 

^y2, Vi) = ^\ (-00, y^; y2±)C{{-oo; fax, yi±)^\\ (-oo, Vi', Via.) (27) 
is the complete gauge link that contains both the transverse [20] 



C-^i-oo] y2±, yiA_) = Pexp 



yi± 

-ig J ■ A±(-oo,£i 



J/2X 



(28) 



and longitudinal gauge link 



>C||(-oo,y ■,y±) = Pexp 



— OO 

-ig J dCA+{C,y±) 



(29) 



The exchange of unphysical gluon between the propagating quark and the 
nucleus therefore leads to the gauge links in the parton matrix elements while 
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physical gluons lead to the higher-twist contributions to the DIS cross section 
which are characterized by the covariant derivatives Dp in the higher-twist 
parton matrix elements. The spatial derivative in Dp comes from the coUinear 
expansion of the partonic hard parts while the transverse gluon field in Dp 
corresponds to interaction between the propagating quark and a physical gluon 
from the nucleus. 

We refer the above organization of semi-inclusive DIS cross section within 
the coUinear expansion as a generalized twist-expansion in which each con- 
tribution involves a parton matrix element with a given number of covariant 
derivatives. These parton matrix elements correspond to unintegrated parton 
distributions. Note that after integration (over both k~ and /c^), each of these 
unintegrated parton distributions will give rise to a mixture of coUinear (or 
integrated) parton distributions or parton matrix elements with different di- 
mensions, because of the off-shellness and transverse momentum carried by 
each parton. 

The nuclear enhanced twist-four contributions to the semi-inclusive cross sec- 
tion at the leading order is contained in Eq. (!23l) with the corresponding 
projected parton matrix elements ujP'uj^'^^^l^,{ki,k2,k) in Eq. fl26|) which is 
apparently gauge invariant. Note that the leading components of the pro- 
jected covariant derivative Dpi are the transverse ones, D_l = + igA^. 
Using the following identity 



idy^C{yi,y) = C{yi,y)[iD±{y ,y±) 



(30) 



it is easy to see that the parton matrix element uj^' uj^' ^^^]^f{ki, ^2, k) in Eq. (12^ 
contains exactly the same quark-gluon correlation distribution in Eqs. (jlj) 
and ([8]). Note that the hard part in Eq. fl23l) is projected to the transverse 
polarization of the gluon exchange and therefore corresponds to interaction 
with coUinear physical gluons in the light-cone gauge. In the derivation [19] of 
this final form. Ward indentities are also used to relate it to the derivative of 
the hard parts for quark and longitudinal gluon interaction in the covariant 
gauge. 

After integration over the transverse momentum of the quark, one will be 
able to recover the nuclear enhanced twist-four contributions to the inclusive 
DIS cross section [22][23] . which should be power-suppressed by 1/Q^. In the 
calculation of the above twist-four contribution to the inclusive DIS cross 
section, one can similarly relate the transverse gauge field A± to the field 
strength F^± via partial integration. However, one should include the gluonic 



13 



poles [16] in 

/ dy-e'^P^y-A^{y-) = [ rfy-e"^^^" (31) 
J J X ±te 



depending on the boundary condition for A±{±oo) which do not vanish si- 
multaneously in the light-cone gauge [23]. We have chosen A±{oo) = in this 
paper. Such gluonic poles have non-vanishing contributions to higher-twist 
inclusive cross sections of DIS which are suppressed by l/Q"^. Without inclu- 
sion of these gluonic pole contributions, one could be misled to the conclusion 
[25] that interaction with physical transverse gluons in the light-cone gauge 
could lead to leading twist transverse momentum broadening of a propagating 
quark, which does not contribute to higher- twist (1/Q^ suppressed) inclusive 
cross section. As we will show in the next section, such transverse momentum 
broadening comes from interaction with unphysical gluons (either longitudinal 
gluon in covariant gauge or pure transverse gauge field A±{—oo) at the 
light-cone infinity in the light-cone gauge). 



5 Nuclear transverse momentum broadening 



Even though the exchange of unphysical gluons do not contribute to the hard 
parts of higher-twist DIS cross sections, they do lead to important gauge 
links in parton distributions and the transverse momentum broadening of the 
propagating quark inside a nucleus at the leading twist (no 1/Q^ suppression) . 
Exchange of physical gluons (e.g. transverse gluons) only leads to higher- 
twist inclusive cross section and transverse momentum broadening which is 
suppressed by factors of l/Q"^ as compared to the leading twist broadening. 
This has been a source of confusion in previous studies. 

Consider the semi-inclusive hadronic tensor with no physical gluon exchange 
in Eq. fl2T]) . The coUinear component of the parton matrix element 



^'^'Hk) = \i>m) + - x^)fUk) (32) 



will give rise to the leading twist contribution to the semi-inclusive hadronic 
tensor 



_ 1 f ,,2, d2)a Z \^^\£i(o)(^\^ 



dH I 2tt 



J dxd'kj^'\l^ - ^^)Tr[^W(x)|]/l(x, k^)- (33) 



and the transverse momentum dependent quark distribution function 
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X {A I V^(0,0^)^£(0,y)^(r,y±) I A). (34) 
The corresponding integrated quark distribution function is 



/l(x)= J d^kUlix^k^) 

= I ^e-^"^-(A I ^(0,0^)^/:||(0,r)V^(l/-,0^) I A). (35) 

One can make a Taylor expansion of the gauge hnk and the quark field in the 
transverse coordinate and then complete integration over y^, 



X 



idv , di- 



-£(0,i/)V;(l/- yl)] \A)5^^\k^). (36) 



Using the identity in Eq. ( !30l) . one can express the above transverse momen- 
tum dependent quark distribution in terms of coUinear higher-twist quark and 
gluon correlation matrix elements |21j . 



For the purpose of discussion in this paper, we consider first the quadratic 
term in the expansion of the gauge link in the covariant gauge, 

£(2)(o,?/) = y dy^dy:^ [-g^A+iy^ ,y±)A+{y^ ,y^)9{y- - y2)9{y2 - Vi) 
+ g^A+iy^, 0±)A+iy^, y±)^(-y2~)^(r - Vi) 
- g'A+{y^, Ox) A+(yr, Ox)^(-yr)^(2/r - 2/2~)] • (37) 

The three terms in the above expansion correspond to the left, central and 
right cut diagrams in Fig.[3]|^c). One can further expand the derivative operator 
exp{idy^ ■ dk^) to the quadratic term. It is easy to note that only the left-cut 
contribution has a term F+_|_(7/^)F+_L(?/f ) ^ f?±^+(l/2^, 0x)(9j_A+(|/f , 0^) from 
the quadratic derivative {dy^ -dk^)"^ ■ As we will explain later this is the leading 
term that has a nuclear enhancement. 

One can also express the above contributions explicitly in terms of the trans- 
verse momentum of each gluonic field as denoted in Fig. E|^c), 



_ ikij_-yi±-ik2±-y2± 

(27r)2 (27r) ' 
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D = J dy^dy^ [-g^e{y~ - y2)e{y2 - yi)5^^\i^ -k^- ki± - hi.) 
+g^{-y:^)e{y' - y^)6^'\i^ - k^ - h±) 

-g^e{-y^)e{y^ - y2)6^'\i^ - k^)] . (38) 

One can expand the above (5-functions in the transverse momenta of the quark 
and gluon fields k±, kij_ and k2±. Contribution to d±A^{y2 ,0±_)d±A+{yi ,0j_) 
again comes only from the first term which corresponds to the left-cut di- 
agram in Fig. [3](c). Since the dependence of the partonic hard part D on 
the transverse momentum of the quark field only leads to the higher-twist 
matrix elements that are not nuclear enhanced, one can set fc^ = in the 
hard part. By momentum conservation ki± = —k2±. In this case, only the 
second term from the central-cut diagram contribute to the parton matrix 
dj_A+{y2 ,0±)dj_A+{yY ,0±) for the leading nuclear broadening. 

One can make a similar analysis of the transverse momentum distribution in 
the light-cone gauge. The leading twist contribution comes from the transverse 
gauge link which is determined by the boundary condition of the transverse 
gauge field A±{—oo,y±) [2T]. As we have discussed before, quark interaction 
with physical transverse gauge field will lead to higher-twist inclusive cross 
section and the transverse momentum broadening which is power suppressed 
by as compared to the leading twist contribution. 

Taking the quadratic derivatives of the two-gluon contribution [Fig. ^c)] to 
the gauge link in Eq. ( l36l) . one can obtain the corresponding transverse mo- 
mentum distribution, 



where 



r,1(^,^i) = / -^dy^dy2e'''^^y-+^^^p^'-y^~y^^e{-y2)e{y^-y^) 

^{A I mi^F+±iy2) ■ F+x(2/r)^(y") I A), (40) 

is the quark-gluon correlation function inside the nucleus. 

The quark-gluon correlation function inside a nucleus in Eq. fl40p has two types 
of contributions. Because a nucleus consists of nucleons which are color singlet 
states, the quark-gluon pair could either come from a single nucleon or from 
two separate nucleons. In the first case, all four parton fields in the above cor- 
relation matrix element are confined to the size of a nucleon y~,yi ,y2 ^ 'T'n- 
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On the other hand, if quark and gluon fields are confined to two separate nu- 
cleons, y~ ,\yi —y2\ ^ ^N, the overall position of the gluon field will follow the 
second nucleon and are only confined to the size of the nucleus Ra- Therefore, 
the quark-gluon correlation function in this case will have a nuclear enhance- 
ment of the order Ra/tn ~ A^^^'^ as compared to the first case where both 
quark and gluon fields are confined to a single nucleon. For multiple scatter- 
ing in a large nucleus, we will only keep the second correlation with nuclear 
enhancement. If we further neglect the correlation between different nucleons 
and assume the large nucleus as a weakly bound and homogenous system of 
nucleons, the leading contribution to the quark-gluon correlation function can 
be approximated as [9)126] 



x(iV|F+x(0)-F+^(r)|iV) 

= 7r/l(x) I d^^pi{^N)xJU^i), (41) 
where pj^ = p'^ is the longitudinal momentum per nucleon and 

xfAx) = J I F+±i^) ■ F+^in I N), (42) 



is the gluon distribution function in a nucleon, respectively. The spatial nu- 
cleon density inside the nucleus is 



The integration over the spatial position of the nucleon = iVi + 2/2 )/2 is 
limited to the size of the nucleus. The relative coordinate of the two gluon 
fields is = - y^ ■ 

With the above approximation, the transverse momentum broadening squared 
of the propagating quark can be calculated as 

(A^l) ^ J d'i^ilnixj^) = J rfG^F(e^), (43) 



where the quark transport parameter 

QFi^N) ^ ^pi{^N)[xf^ {x)Uo (44) 
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can be interpreted as the broadening of the mean transverse momentum 
squared per unit path length. One can consider multiple gluon exchanges the 
transverse distribution will become a Gaussian form pT]|25f29] with the width 
given by the averaged transverse momentum broadening. 

In the above approximation of the twist-four quark-gluon matrix we have ne- 
glected multiple-nucleon correlation in a large nucleus. Such an approximation 
is not valid for small x where quark-gluon and gluon-gluon fusion from dif- 
ferent nucleons become important and can lead to modification of the quark 
distribution function and gluon saturation in a large nucleus P7f25] . One can 
take into account such effect by using a nuclear modified quark distribution 
function f\{xB) and saturated gluon distribution function in the transport 
parameter qp which could lead a non-trivial nuclear and energy dependence 

Note that the quark-gluon correlation function as defined in Eq. fl4Up and the 
gluon distribution function in Eq. fl42l) are not gauge invariant in covariant 
gauge. One needs to resum additional number of coUinear soft gluons on both 
side of the cut to produce gauge links that will ensure the gauge invariance 
of the quark-gluon correlation and the gluon distribution function. A general 
and gauge invariant form of the transverse momentum broadening has been 
derived in Ref. [2T] . 



6 Induced gluon radiation in light-cone gauge 

As another example of the equivalence of the hard parts in double parton 
scattering in light-cone and covariant gauge, we consider the double hard 
quark-gluon scattering in next leading order which has been calculated in co- 
variant gauge in helicity amplitude approximation in Ref. [8] . In the following 
we will calculate the induced gluon spectra from such quark-gluon scattering 
in light-cone gauge approach within helicity amplitude approximation. Under 
such helicity amplitude approximation, we can neglect momentum transfer to 
the quark, except in its propagation direction and only consider its dominant 
minus component. Since the initial gluon fields only have transverse compo- 
nents in the light-cone gauge approach, their direct interaction with the quark 
will produce a vertex in the form r/i^±_i/i = in the helicity amplitude approx- 
imation. Intuitively, this is because a quark can not absorb the transversely 
polarized gluon due to helicity conservation if we neglect the recoil induced 
by the interaction. For this reason, we only need to consider the quark-gluon 
rescattering with triple gluon vertex in the helicity amplitude approximation 
as shown in Fig. HI Inclusion of the quark recoil and other diagrams in light- 
cone gauge will lead to power corrections in i.\ and 1 — z. 
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Fig. 4. Radiative correction to double hard scattering 

Contribution from the next-leading order double hard quark-gluon scattering 
can be written as, 

dz gJ^TT J 271 2tx 

x/^H^;::{x,xuX2,p,q,z){A \ ^,(0)^F+.(2/2-)F/(yr)V^,(y") I ^) ,(45) 

where X3 = x + Xi — X2 by momentum conservation and z is the fractional 
momentum carried by the final quark i~ = zq~ . The hard partonic part in 
general has the form, 



HP:ix,x^,X2,p,q,z)= J -0-^27r5^{f)6{l-z-^—) 

x^Tt[^j,HP-^,]. (46) 
For the central-cut diagram of quark-gluon rescattering in Fig. HI H^"^ is given 

by, 



{xp + qy — ie {x^p + qY + ie {xip — £y — ie 



^((^T^T^^^F^T)^^""'^^^'^'^'^ ^ 

{2p ■ qY{x — xb — — xb + i^) (1 — zY 

[xi + ie)[x + Xi — x-^ — ze) 
where C = Ca/2Nc is the color factor, £ and iq = {xi + X2)p + q — £ are the 
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momenta of the final state gluon and quark, respectively, and 



is the fractional momentum taken away from the initial state quarks and 
gluons by the final state quark-gluon splitting. The three gluon vertices are 
defined as 

j^a'a^ = ^"'"(2x2P - iy - g'^^ii + X2vT' + ^^"'(2£ - x^vY ■ (49) 

In the above, we have used gluon propagators in the Feynman gauge. Since 
the initial gluons only have transverse components in their polarization, we 
can also replace the summation of the final state gluon's polarization tensor 

One can carry out the integrations over x and x\ in Eqs. (H5l) and (H7|) by con- 
tour integration. There are four possible poles in the denominator of Eq. (HTIl 
from quark and gluon propagators. Different choices for the pair of poles rep- 
resent subprocesses with different kinematics. Here we choose the poles at 
X = xb,Xs = xb which correspond to the double hard scattering. In this case 
the momentum fraction carried by the initial gluon is Xi = xl which scatters 
with the propagating quark that has initial momentum xb'P + = [0, g~, 0_l]. 
Therefore, in this case of double hard scattering, the hard part can be written 
as 



HP- = (xB^ + ^)_LTr[7-if^-] . (50) 

Since the contributions from soft gluon radiation {z 1) is dominant in 
fragmentation process, we can take the helicity amplitude approximation, /g ^ 
zq~r/t in the trace of the hard part. We have then. 



dx 



dxi dx2 
4 



d. 



pa 



2zq- 



2xiX2 

{xBi) + ^)27r2j9 ■ q 



HP"" 



{2p ■ g)4(l - ^)V 



Tr 



d, 



epg' 



-[{xB^ + i)27i2p-q]. 



(51) 



For the soft radiation approximation z ^ 1, one obtains in this light-cone 
gauge approach. 
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^^^i:^J>.PW)|^^T,-,(.....). (52) 

which is the same as the resuh for double hard scattering in the covariant 
gauge [8]. 

In the calculation of semi-inclusive DIS, the transverse momentum enters 
as another scale in additional to of the virtual photon. The twist-four 
contribution to the semi-inclusive spectra at the order of result of 

the induced gluon bremsstrahlung is suppressed by as/(-\ as compared to the 
leading twist result, 



The coUinear expansion in this kind of inclusive processes is only valid when l\ 
is much larger than the average transverse momentum of the initial parton or 
the quark transverse momentum broadening which is given by the jet transport 
parameter q [Eq. fHSl) ]. We have also neglected contributions proportional 
to for l\ <^ Q"^ . For small values of l\ <^ qRA the twist-expansion 

method will fail for the semi-inclusive processes and one needs to regularize 
the divergency of the semi-inclusive spectra. However, the LPM interference 
between double hard and soft rescattering processes suppresses the induced 
spectra for small £^i?^/2g^ <^ 1. Therefore, the final result in the coUinear 
expansion will be a good approximation and insensitive to the regularization 
for large initial quark energy q~ ^ qR\- 



7 Summary 



In this paper, we have investigated the gauge invariance of the leading twist- 
four contribution to the semi-inclusive cross section of DIS off a large nucleus 
due to multiple parton scattering in the framework of generalized coUinear 
factorization. 

We first proved the general equivalence of the hard parts of double scattering 
in light-cone and covariant gauge, using a set of identities for hard partonic 
processes which were derived from Ward identity and equation of motion. 
This equivalence hold for double parton scattering in the twist expansion to 
all orders in «s. We also give two specific examples to demonstrate the equiv- 
alence explicitly. It is easy to see that our proof can be directly extended to 
other higher twist contributions. The equivalence of the hard parts in different 
gauge at any twist level can be proved by connecting the contributions of the 
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transverse and longitudinal components of the gluon using Ward identity as 
we did. 

We have also demonstrated explicitly the gauge invariance of higher twist 
contributions in the calculation of semi- inclusive DIS cross section in the lowest 
order and the next-leading-order with induced gluon emission. We pointed 
out the importance of the gluonic poles in the calculation of the higher-twist 
contributions and that interaction with transverse gluons only leads to higher- 
twist inclusive cross sections which are power suppressed by l/Q^. The leading 
contribution to transverse momentum broadening comes from interaction with 
only the unphysical gluons. 
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